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Abstract 

Analytic-bilinear approach for construction and study of integrable 
hierarchies, in particular, the KP hierarchy is discussed. It is based 
on the generalized Hirota identity. This approach allows to represent 
generalized hierarchies of integrable equations in a condensed form of 
finite functional equations. Resolution of these functional equations 
leads to the r-function and addition formulae to it. General discrete 
transformations of the r-function are presented in the determinant 
form. Closed one-form and other formulae also arise naturally within 
the approach proposed. Generalized KP hierarchy written in terms of 
different invariants of Combescurc symmetry transformations coincides 
with the usual KP hierarchy and the mKP hierarchy. 

1 Introduction 



The Kadomtsev-Petviashvili (KP) equation and the whole KP hierarchy of 
equations are significant parts of the theory of integrable equations. They 
arise in various fields of physics from hydrodynamics to string theory. They 
are also the tools to solve several problems in mathematics from the differ- 
ential geometry of surfaces to an algebraic geometry. 
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The KP hierarchy has been described and studied within the framework 
of different approaches. The Sato approach j| (see also @-@) and the 8- 
dressing method ||-|]|] are, perhaps, the most beautiful and powerful among 
them. The infinite-dimensional Grassmannian, pseudo-differential opera- 
tors, Hirota bilinear identity and r- function are the basic ingredients of the 
Sato approach which is in essence algebraic. In contrast, the <9-method, 
based on the nonlocal 9-problem for the wave function, mostly uses the 
analytic properties of the wave function. These two approaches look like 
completely different. On the other hand, each of them has its own advan- 
tages. So one could expect that their marriage could be rather profitable. 

A bridge between the Sato approach and the 9-dressing method has 
been established by the observation that the Hirota bilinear identity can 
be derived from the 5-dressing method |§, Elements of the approach 
which combines the characteristic features of both methods, namely, the 
Hirota bilinear identity from the Sato approach and the analytic properties 
of solutions from the 9-dressing method have been considered in [1C], [11 1. 

This paper is devoted to the analytic-bilinear approach to integrable hi- 
erarchies. It is based on the generalized Hirota bilinear identity for the wave 
function with simple analytic properties (Cauchy-Baker-Akhiezer (CBA) 
function). This approach allows us to derive generalized hierarchies of in- 
tegrable equations in terms of the CBA function ip(\, fj,, x) in a compact, 
finite form. Such a generalized equations contain integrable equations in 
their usual form, their modified partners and corresponding linear problems. 

The generalized Hirota bilinear identity provides various functional equa- 
tions for the CBA function. The resolution of some of them leads to the 
introduction of the r-function, while others give rise to the addition theorems 
for the r-function. The properties of the r-function, such as the associated 
closed 1-form and its global definition are also arise in a simple manner 
within the approach under consideration. Generic discrete transformations 
of the CBA function and of the r- function are presented in the determinant 
form. These transformations are the generalized form of the usual iterated 
Darboux transformations. 

The generalized KP hierarchy is presented also in the 'moving' frame 
depending on the parameter. This generalized KP hierarchy written in the 
different 'moving' frames contains the Darboux system of equations. 

In addition to the usual infinite-dimensional symmetries the generalized 
KP hierarchy possess the symmetries given by the Combescure transforma- 
tions. The invariants of these symmetry transformations have the compact 
forms in terms of the CBA function. The generalized KP hierarchy written 
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in terms of these invariants coincide with the usual KP hierarchy , mKP 
hierarchy and SM-KP hierarchy. 

The present paper is devoted to the one-component KP hierarchy. The 
authors plan to consider multi-component KP hierarchy and the 2-dimensional 
Toda lattice in subsequent papers. 

The paper is organized as follows. Generalized Hirota identity is intro- 
duced in section 2. The generalized KP hierarchy is derived in section 3. 
Combescure symmetry transformations are discussed in section 4. Gener- 
alized KP hierarchy in the 'moving' frame is considered in section 5. The 
r-function is introduced in section 6. The addition formulae for the r- 
function are also obtained here. Transformations of the CBA function and 
the r-function given by the determinant formulae and the Darboux trans- 
formations are discussed in section 7. Closed one-form variational formulae 
for the r-function are presented in section 8. 

2 Generalized Hirota identity 

The famous Hirota bilinear identity provides us condensed and compact 
form of the integrable hierarchies. Here we will derive the generalized Hirota 
bilinear identity in frame of the 9-dressing method. 

The 9-dressing method (sec || - ||) is based on the nonlocal <9-problem 
of the form 

a A ( X (x, A) - 77(x, A)) = JJ dfiA dfix(M)9~Hl*)Rfa U) 

(x(x, A) - r?(x, A))| AHoo -> 0. 

where A € C, d\ = d/dX, r/(x, A) is a rational function of A (normalization). 
In general case the function x(^) an d the kernel R(X, //) are matrix-valued 
functions. 

The dependence of the solution x(A) of the problem (Q) on the dynamical 
variables is hidden in the function g(X) . Here we will consider only the case of 
continuous variables, for which gi = exp(JQxj), where Ki(X) are, in general, 
matrix meromorphic functions. 

It is assumed that the problem (|]) is uniquely solvable. The 5-dressing 
method allows us to construct and solve wide classes of nonlinear PDEs 
which correspond to the different choice of the functions Ki(X) . 

Here we will assume that the kernel R(X, fx) is equal to zero in some open 
subset G of the complex plane with respect to A and to ji. This subset should 
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typically include all zeroes and poles of the considered class of functions g(X) 
and a neighborhood of infinity. 

In this case the solution of the problem (1) normalized by 77 is the function 

X(A) = r/(x, A) + <p(x,A), 

where 77(A) is a rational function of A (normalization), all poles of 77(A) 
belong to G, tp(X) decreases as A — > 00 and is analytic in G. 

The special class of solutions of the 3 problem (||) normalized by (A— ^)~* 
(77(A) = (A — ju) _1 , fj, G G) is of particular importance for the whole <9- 
dressing method ||. Let us consider the 9-problem (|I|) for such solutions 
and corresponding dual problem for the dual function x*(A,/u): 

9ax(A,ju) = 2m5(X - fi) + J J du A dvxiy^)gi{v)R{v, \)gi{uy l , 
a AX *(A,M) = -27ri5(A-M)- JJ c dvAdvg 2 {v)R{\v)g 2 {v)- l X *{v,n). (2) 
where #i = <7(A,x), 52 = g(X,x') After simple calculations, one gets 
JJ^duAdu-^ (x{v, \; gi)g 1 (v)g 2 (v)~ 1 x* (v, V, g2)) = 

/ xiy, ^,gi)gi(v)g2(vy 1 x*(v,Li;g2)di< = o. (3) 

In the particular case g\ = g 2 from (||) it follows that in G the function 
x(A, //) is equal to — x*(/x, A) (see also li~2(] ). So finally we have 

/ x(^WS'i)5i(%2 _1 ( l/ )x(A, v;g 2 )dv = 0. (4) 
•/<9G 

Here the function x(A,/j) possesses the following analytical properties: 

<9ax(A,aO = 2vri(5(A - /i), -<9 M x(A,/i) = 2vriJ(A - /i), 

where <5(A — fi) is a (^-function, or, in other words, x - > (A — A*) -1 as A — > /i 
and x(A, A 4 ) is analytic function of both variables A, \i for X ^ fi. 

In the particular case A = fi = the relation (||) is nothing but the usual 
Hirota identity for the KP wave function x{ v i 0, x) and the dual KP wave 
function %*(z/, 0,x') (see e.g. |]]-[§|])- Just in this case (for A = \i = co ) the 
Hirota bilinear identity has been derived from the 3-problem in |J. 

The identity @ represents itself the generalization of the Hirota identity 
which is bilocal with respect to the dynamical variables x and the spectral 
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variables (A,/i). It can be considered as the point of departure without any 
reference to the <9-dressing method. Namely, the generalized Hirota bilinear 
relation (||) is the starting point of the analytic-bilinear approach which we 
will consider in this paper. The double bilocality (with respect to x and 
(A, fi) provides us an additional freedom which will allow us to represent 
integrable hierarchies in a unified and condensed form. 
Introducing the function tp(X,fi,x) via 

VK^ /■*;#) = g~ l {p)x{\vig)gW , (5) 

one gets another form of the generalized Hirota equation 

^{v^]gi)^{\v;g 2 )dv = 0. (6) 

dG 

Note that in the framework of algebro-geometric technique the function 
V>(A, jj) corresponds to the Cauchy-Baker-Akhiezer kernel on the Riemann 



surface (see [ 13 1 ) . We will refer to the function fi) as the Cauchy-Baker- 
Akhiezer (CBA) function. 

We will assume in what follows that we are able to find solutions for 
the relation (Q) somehow. In particular, it can be done by the 9-dressing 
method. 

The general setting of the problem of solving (||) requires some modifica- 
tion of Segal- Wilson Grassmannian approach Q . Let us consider two linear 
spaces W(g) and W(g) defined by the function x(A,/i) (satisfying (Q)) via 
equations connected with equation @ 

f(v;g)x{\v;g)dv = 0, (7) 

dG 

X{u,n;g)h(u;g)du = 0, (8) 

dG 

here /(A) G W, h(X) G W; /(A), h{\) are defined in G. 

It follows from the definition of linear spaces W, W that 

/(A) = 2vri jj r,(u) X (X, v)dv A dP, n(u) = (J=f^)) , 

h(ji) =- 2mjj G x(v,»Wv)dvAdv, rj{v) = (SzK^i) ■ ( 9 ) 



These formulae in some sense provide an expansion of the functions / , h 
in terms of the basic function x(A, fj)- The formulae (0) readily imply that 
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linear spaces W, W are transversal to the space of holomorphic functions in 
G ( trans versality property). 

From the other point of view, these formulae define a map of the space 
of functions (distributions) on G r], rj to the spaces W, W. We will call rj 
(rj) a normalization of the corresponding function belonging to W (W). 

The dynamics of linear spaces W, W looks very simple 

W(g)=W g~ 1 ; W(g) = gW , (10) 

here Wq = W{g = 1), Wo = W(g = 1) (the formulae 
identity (Q) and the formulae @). 

To introduce a dependence on several variables (may be of different type), 
one should consider a product of corresponding functions g(\) (all of them 
commute). 

The formula (01) is a basic tool for our construction. Analytic properties 
of the CBA kernel accompanied with the different choices of the functions 
<7i and gi will provide us various compact and useful relations. 



10) follow from 



3 Generalized KP hierarchy 

In the present paper we will consider the scalar KP hierarchy. It is generated 
by the generalized Hirota formula (|j) where G is a unit disk and 

5 (x,A)=exp(f>A-M • (11) 



\i=l 



Let us consider the formula (0) with 



9i9 2 1 = 5(x, v)g ^x', v) = exp \ Y^(xi ~ x\)v l j = ^— ^ . (12) 



where a and b are arbitrary complex parameters, a ,b G G. Since log(l — e) 
Ya^=i e *A> one has 

1 1 

x'j — Xi = -a 1 b l . 

i i 

Substituting the expression ( |l2"| ) into (||), one gets 

; ' ^ X(A, x + [a]) - (\ — ^] x(A, pi, x + [b]) + 



fx — b J \X — b , 

(6-a)x(A,6,x+[a])x(6,M,x+[&]) = 0, A^M (13) 
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where x+[a] = Xj + [a]j ,0 < i < oo , [a]i = \a l . Equation (13) is the simplest 
functional equation for %(A,//,x) which follows from the generalized Hirota 
equation (f|). 

Residues of the l.h.s. of (|i~3|) at the poles /i = b and A = b vanish 



identically. Evaluating the l.h.s. of (12) at fi = a and A = a, one gets the 
equations 



J x(A, a, x + [b]) = (b- a) X (X, b,x + [a}) X (b, a, x + [b]) , (14) 
X(a, //, x + [a]) = (a- b) X {a, b, x + [a]) X {b, /x, x + [b]) , (15) 



A — a 
X-b 
fi — a 
fx — b 

a^b. 

Equations ([l4|) and (|l5|) imply that 
(A - /z)x(A,/i,x) = 

(A - a)x(A,q,x - [a] + [/i]) = (a - ^)x(a,^,x + [a] - [A]) 
(A* - °M/A a, x - [a] + [//]) (a - X)x(a, A, x + [a] - [A]) ' 
Since (/i — a) X (fj,, a) — ► 1 as /i — a — > 0, one gets from (p^) 

(A - /i) 2 X (A, /x, x + [X]) X (ji, A, x + [/i]) = -1 . (17) 



We will solve the equations ( |14|) -(|17|) in the next section. Now, let us 
consider the particular form of equation (|l3|) for 6 = 0. In terms of the CBA 
function it reads 

V?(A, /j,x+[a])— V>(A,it,x) = a-0(A, 0,x+[a]W0,//,x); x'-Xi = ^a*. (18) 

This equation is a condensed finite form of the whole KP-mKP hierarchy. 
Indeed, the expansion of this relation over a generates the KP-mKP hierar- 
chies (and dual hierarchies) and linear problems for them. To demonstrate 
this, let us take the first three equations given by the expansion of ( |T8| ) 
over a 

a: ip(X,fj,,x) x = ^(A,0,x)V>(0,//,x) , (19) 

a 2 : V(A, /J, x) y = ip{X, 0, x) x tp(0, fi, x) - V(A, 0, x)^(0, fi, x) x , (20) 

1 3 

ip(X, n,x.) t = ~tp(X,fj,,x) xxx - -ip(X,0,x) x -(p(0, n,x) x + 

j (V(A,0,x)^(0,/x,x) - ^(A,0,x)^(0,m,x) v ) (21) 
x = xi; y = x 2 ; t = x 3 . 



,3 
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In the order a 2 equation (^) gives rise equivalently to the equations 

ip(X,fj,,x) y -ip(X,fi,x) xx = -2^(A,0,x)V>(0, fj,,x) x , (22) 
ip(X,(i,x)y + ip(X, p,x) xx = 2-0(A,O,x) :r V(O,/U,x) , (23) 

Evaluating the first equation at p = 0, the second at A = one gets 

/(x)„ - fjx) xx = u(x)/(x), (24) 
7(x)„ + /(x)^ = -n(x)/(x) (25) 

where u(x) = —2-0(0, 0) x and 



f = J ^(X,0)p(X)dX. 



p(X) and p(p) are some arbitrary functions. 

In a similar manner, one obtains from (|l~9|)- (^j) the equations 

3 3 

/t - fxxx = -juf x + -(u x + d~ u y )f , (26) 
~ ~ 3 ~ 3 ~ 

ft ~ fxxx = -xuf x + -(u x - d~ 1 u y )f . (27) 



Both the linear system (24), (pll) for the wave function / and the linear 



system (|25D, (|27D for the wave function / give rise to the same KP equation 

I 3 3 

~^U XXX + — UU X T ~J 



a i ~ -Uxxx + ~uu x + -d x 1 u yy . (28) 



To derive linear problems for the mKP and dual mKP equations, we inte- 
grate equations (|l~9|) , (|22]), (|23|) and (^T|) with the two arbitrary functions 
P(A), p(aO 

^(x), = /(x)/(x), (29) 
<&(x) y - $(x) OT = -2/(x)/(x) x , (30) 
$(x) y + $(x) aa: = 2/(x) x /(x), (31) 

SCx^-SCx)^ = -|/(x) x /(x) !B -|(/(x)7(x) 1/ -/(x) 1 ,/(x))(32) 
where 

$= // p(n)i>{X,n)p(X)dXdn. 
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Using the first equation to exclude / from the second (and / from the third) , 
we obtain 

$y-<$>xx = V(x)$x, (33) 

$ w + *xx = -v(x)* a (34) 
where v = -24£k, v = 2^Pf. 



Similarly, one gets from ( |2l| ) 

3 3 
= -^(x)$ xx + -(v x +v 2 + d-\)$ x , (35) 



3 3 

$t - = ^W*** + 4^ + y2 ~ 3* %)** • (36) 



The system (|S3), (S3) gives rise to the mKP equation 



3 

''/ = Vxxx + J v2v x + %v x d~ v y + 3d~ v yy , (37) 



while the system (|34), (|36 ) leads to the dual mKP equation, which is ob- 
tained from the ( |37[ ) by the substitution i> — > v, t — > — t, y — > — y, x — > — x. 

So the function is simultaneously a wave function for the mKP and 
dual mKP linear problems with different potentials, defined by the dual KP 
(KP) wave functions. 

Using the equation ( |2l| ) and relations ( |33|) and ( |34[ ) , one also obtains an 
equation for the function <I> 

1 3 <3? 2 — <3? 2 3 <]? 

$t "f«"8^r + 4^ = °' ^ = i;- (38) 

This equation first arose in Painleve analysis of the KP equation as a sin- 
gularity manifold equation ||14[| . 

It is tedious but absolutely straightforward check that the expansion of 
( |l~8| ) in higher orders of a generates 

1) the whole hierarchy of KP singularity manifold equations for ^(A,//) (or 
<Kx) 

2) the hierarchy of linear problems for the mKP and dual mKP equations, 
where ijj(\,fi) (or$(x) is the common wave function and v = — 2(log^(0, A, x))^, 
v = — 2 (log ip(X, 0, x))a; are the potentials 

3) mKp hierarchy for v and dual mKP hierarchy for v 

4) the hierarchies of KP linear problems for ?/>(A,0,x) and dual KP linear 



9 



problems for ip(\,0,x) 
and, finally 

5) the KP hierarchy of equations for u = —2vp(0,0) x . 

Note also one interesting consequence of the formula (|j 

x(0,A,x) = -^x _1 (A,0,x + [A]). 



(39) 



4 KP hierarchy in the 'moving frame'. Darboux 
equations as the horizontal subhierarchy 



Now let us consider the expansion of the l.h.s. of fll3|) over e = a — b, where 
e — > 0. In the first order in e one gets 



where 



Ai(6)V>(A,/i,x) = ^(&,/z,x)V>(A,&,x) 

un-l 9 



(40) 



A 1 (b) = E 



dx r . 



n=l " J " n 

In the higher orders in e one obtains the hierarchy of equations of the 
form (|i~9D-(f2i~D and their higher analogues with the substitution ^(A, 0, x) — > 
ip(X,b,x), ip(0, /x,x) — > ip(b, fx, x) and 



_d_ 

dxi 



Mb) = E 



nl 



'. n(n — dxi 



Such a substitution is in fact nothing but the change of dynamical variables 
(or the coordinates on the group of functions g). Indeed, it is not difficult 



to show that Aj(&) 
by the relation 



It is clear that 



aa 5 , where the dynamical variables Xi(b) are defined 
Xi{b) 



E 



E 



d d 



dxiix'Y dxi(\) 



0. 



Note one interesting property of the derivatives g^ro , namely 



5 a 



<9A' dxi(\) 



(i + l) 







dxi+i(X) 



(41) 



(42) 



(43) 
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So the operator ^ is a 'master symmetry' for all vector fields q^J)- 

The expansion of equation ( |l3| ) up to the third order in e gives the 
equations 

Mi x(b)) = i>(b, n, x(6))V(A, 6, (x(6)), (44) 



dx2{b) dxi(b) dx\{b) 

d Id 3 3d d 

-1p(X, H, X(b)) = - 77^- 0(A, M) - T o /,. ^(A, b) ■ j7T ip(b, n) + 



dx 3 (6) rv x " Ad Xl (bY YV ^' 4dxi(by v ' J d Xl {b) 

+ 1 ,MK 61 ' * ( "' M) " **• « ■ a£® M) ) • (46) 

The analogues of equations (|2^), (^) have the form 



-^(A, M ,s(6))-— -7- -^(A, M )+2^(A,6)— — -^(6, M )=0, (47) 



dx2(b) ' ' Oxi(6) 2 ' ' dx\(b) 



' ^{\, l i,x{b)) + — ip(\, fI )-2— — ^(\,b)-^(b,(J,)=Q. (48) 



dx 2 (b) Yy wy 9xi(6) 2rv 8 Xl (b) 
Equations (44)- (46) and higher equations again give rise to the general- 



ized KP hierarchy but now in coordinates Xi(b),i = 1,2,3.... For such KP 
hierarchy written in the 'moving' frame the parameter b is an arbitrary one, 
but fixed. 

Let us consider now equations of the type ( |44l ) written for several values 
of b. We denote xi(b a ) = £ Q ,a = 1,2, ...,n. Equations ( f44|) taken for b = 
b a , A = b/3, \i = 6 7 (a ^ (3 ^ 7), look like 

d 

7^t-# 7 = ^/3o'0o7 1 « / + 7 (49) 

OC,a 

where ifi a p = ip(b a ,bpx). The system (E9) is just well-known system of 
N 2 — iV resonantly interacting waves. 

Integrating equations ( fi"4|) over // with the function p(//) and evaluating 
the result at b = b a , 7 = bp, one gets 

^=VW«, («//?) (50) 
Ota 
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where fp = J dpip(b/3 ) p)p(p). Analogously one gets 



df* 

q£ = MZ, (a ±&) (51) 

where fg = J d\ip{\,bp)p*(X) and p*(X) is an arbitrary function. The sys- 
tems (|50|) and (51) are the linear problem and dual linear problem for equa- 



tions fl49D , respectively. 

Expressing ip a p via f a and /*, one gets from (pC|) and ( |5l|) (using (p9|)) 
the same system for f a and /* 

d 2 H a 1 dHfsdH a \ dH 1 dB (X fK0 , 



The system (^) is the Darboux system which was introduced for the first 
time in the differential geometry of surfaces [14] and then was rediscovered in 
the matrix form within the d— dressing method in the paper [5]. Note that 
the Darboux equations in the variables of the type x\(b a ) have appeared 
also in the paper [16] within completely different approach. 



One can treat the Darboux equations (|52j) with different n as the hori- 
zontal subhierarchy of the whole generalized KP hierarchy. Note that equa- 
tions (^4|)-(|46|) and their higher analogues give rise to the higher resonantly 
interacting waves equations. 

5 Combescure symmetry transformations for the 
generalized KP hierarchy 

Let us consider now the symmetries of the equations derived above. All 
the higher equations of the hierarchy £1X6 , £IS usual, the symmetries of each 
member of the hierarchy. Here we will discuss another type of symmetries. 



Since p(X) and p(p) are arbitrary functions, equation (|38|) and the hier- 
archy ( |i"8| ) possess the symmetry transformation 

$(p(\),p(p))^$' = $(p'(\)J(p)). 

This transformation is, in fact, the transformation which changes the nor- 
malization of the wave functions. The fact that such transformations are 
connected with the so-called Combescure transformations, known for a long 
time in differential geometry, was pointed out in pi 
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The Combescure transformation was introduced last century within the 



study of the transformation properties of surfaces (see e.g. IIJ, |i7|| ). It is 
a transformation of surface such that the tangent vector at a given point of 
the surface remains parallel. The Combescure transformation is essentially 
different from the well-known Backlund and Darboux transformations. The 
Combescure transformation plays an important role in the theory of the 
systems of hydro dynamical type pq|. It is also of great interest for the 



theory of (2+l)-dimensional integrable systems [19|. Combescure symmetry 
transformations are essential part of the analytic-bilinear approach. 

The Combescure transformation can be characterized in terms of the cor- 
responding invariants. The simplest of these invariants for the mKP equa- 
tion is just the potential of the KP equation L-operator expressed through 
the wave function 

u = /(x)y ~ {^ )xx , (53) 
/(*) 

u = /(x)j C f ^ x)xx , (54) 
/( x ) 

or, in terms of the solution for the mKP (dual mKP) equation 



v' y + v' xx - -((v') 2 ) x = v y + v xx - -{v 2 ) x , (55) 

v y ~v' xx - -{{v'f) x = v y -v xx --{v 2 ) x . (56) 

The solutions of the mKP equations are transformed only by a subgroup of 
the Combescure symmetry group corresponding to the change of the weight 
function p(fi) (left subgroup) and they are invariant under the action of the 
subgroup corresponding to p(X) (vice versa for the dual mKP). 

All the hierarchy of the Combescure transformation invariants is given 
by the expansion over e near the point x of the relation ( |l8| ) rewritten in 
the form 

d //(x) - /(x')\ 13„, , 1 , 

The expansion of the left part of these relations gives the Combescure trans- 
formation invariants in terms of the wave functions /, /. To express them 
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in terms of mKP equation (dual mKP equation) solution, one should use 
the formulae 

v = -2 f -j, / = exp(-^-y); (59) 

v = 2 f -j, f = exp^d~ l v). (60) 

It is also possible to consider special Combescure transformations keeping 
invariant the KP equation (dual KP equation) wave functions (i.e. solutions 
for the dual mKP (mKP) equations). The first invariants of this type are 

^■(x) <Mx) 

/'(*) /(x) ' 

^■(X) <^(x) 

/'(*) /(x) • 

All the hierarchy of the invariants of this type is generated by the expansion 
of the left part of the following relations over e 



(61) 
(62) 



( * ( *}-* (x) ) = ./(*'), 4-« = Y< (63) 

( $(X) / ^ (X ' ) ) = ,/V), xj-x, = -^. (64) 
Now let us consider the equation (|38|) and all the hierarchy given by 



the relation (18). This equation admits the Combescure group of symmetry 
transformations &(p(\),p(fi)) —><&' = $(p'(A), j?{p)) consisting of two sub- 
groups (right and left Combescure transformations). These subgroups have 
the following invariants 

v = ~ *** (65) 

and 

g= + (66) 



From (|33|), (34) it follows that they just obey the mKP and dual mKP 
equation respectively. The invariant for the full Combescure transformation 
can be obtained by the substitution of the expression for v via ( |65|) to the 
formula (55). It reads 

/cE> \ (b cf> 2 — d> 2 

u = 9 x \-^r) — ^ — I • (67) 
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From (ID, (§3|), (g), (H), (||) it follows that u solves the KP equa- 
tion. 

So there is an interesting connection between equation (|38[), mKP-dual 
mKP equations and KP equation. Equation ( |38|) is the unifying equation. 
It possesses a Combescure symmetry transformations group. After the fac- 
torization of equation (|38| ) with respect to one of the subgroups (right or 
left), one gets the mKP or dual mKP equation in terms of the invariants for 
the subgroup ( p5| ) , (Q) . The factorization of equation (|3^) with respect to 
the full Combescure transformations group gives rise to the KP equation in 
terms of the invariant of group (|67|). 

In other words, the invariant of equation (]38| ) under the full Combescure 
group is described by the KP equation, while the invariants under the action 
of its right and left subgroups are described by the mKP or dual mKP 
equations. 

6 r-function and addition formulae 

Now we will analyze the functional equations (|l3|) -(|r7|). Equation (|l6|), 
evaluated at /i = for some a = clq gives 

w\ n \ (A-a)x(A,a ,x- [q ]) Z(X,x) 

Ax(A, 0, x = — — p-r— = — — — - 68 

{-a)x{v,a,x- [a ]) Z{0,x) 

where we denote Z(X,x) = (A — ao)x(A, ao,x — [ao])- Substituting the 
expression (pq) into equation (|T6|), we get 

It is easy to check that in virtue of (|69| ) equation ( |T^) is satisfied identically, 
while equation (|15|) takes the form 



i?(a, A)i?(A, 6)i2(6, a) = i?(a, &)i?(6, A)i?(A, a) (70) 

where i?(a, 6,x) = Z(a, x + [a] + [6]). In terms of R(a, b) we have 

w > \ R(X,ti,x- [A]) 

A - MJX(A, x = — -— . 71 

R{fj,,fi,x- [fj]) 

Thus the problem of resolving equations (0), ( |l5| ) is reduced to the single 
functional equation (170) , which is of the form of the triangle (Yang-Baxter) 
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equation, well-known in the quantum theory of solvable models (see e.g. 

13). 

Prom the definition of R(a,b,x) it follows that it has a certain special 
structure. Indeed, since Z(a, x) = R(a,b,x — [a] — [b]), one has R(a,b,x. — 
[a] — [b]) = R(a,0,x — [a]). Consequently i?(a, 6, x) = R(a,0,x + [b]). 

So we should solve the triangle equation ( |70[) within the class of R of 
the form i?(a, 6, x) = H a (x+ [6]), where H a (x) is some function. Taking the 
logarithm of both parts of (|70|) , one gets 



0(a, A) + 6(A, b) + 0(6, a) = 0(a, b) + 0(6, A) + 0(A, a) (72) 

where = logi?. Representing G as 0(a, 6, x) = + (a,6, x) + 0_(a, 6, x), 
where 0+ and 0_ are respectively symmetric and antisymmetric parts of 
0, one easily concludes that 0+ solves ( |72[ ) identically while 0_ satisfies 
the equation 

0_(a,A) + 0_(A,6) + 0_(6,a) =0. (73) 
Taking equation ( |73| ) at 6 = 0, one gets 

9_(a,A) = 9_(0,A) -0_(O,a). (74) 

Then since 0(a,6,x) (as i?(a,6,x)) has the form 0(a,6,x) = Z a (x + [6]), 
where Z a are some functions, it follows from (f74|) that 



e_(a,A) =Z _(x+[A])-Z _(x+[a]). 

Then for the symmetric part of + one has Z a+ (x + [6]) = ^t+(x + [a]). 
Taking 6 = 0, one concludes that Z a+ (x) = Zq + (x + [a]) So 0+(a, 6) = 
Zq + (x + [a] + [6] ) Thus general solution of (|72|) has the form 

0(a, 6,x) = Z 0+ (x + [a] + [6]) + Z _(x + [6]) - Z „(x + [a]) . 

Consequently, the general solution of (f70|) for our class of R(a,b,x) reads 



i?(a, 6, x) = i? s (x + [a] + [6]) r(x + W) (75) 

r(x+ [aj) 

where i? s and r are arbitrary functions. Substituting now the expression 
(75) into the expression (71), we get 

1 t(x- IX] + \u]) 

X A, /i, x = Li- m > 76 

(A - /i) r(x) 
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This formula coincides with the formula introduced in the paper [11] in 
a completely different context. Note that in our approach the function 
r(x) is still an arbitrary function giving a general solution of the functional 



equations (14), (15) through the formula ([ 

Now we will use the general equation (|i~3|). Substituting (ff6|) into dl3|), 
one gets 



(a - M )(A - 6)t(x + [a] + M)r(x + [A] + [b]) + 
( A _ a)(b _ M ) r(x + [ A ] + [a ]) T ( x + [6 ] + M) + 

(6 - a) (A - m)t(x + [6] + [a])r(x + [A] + [fx]) = . 



(77) 



It is nothing but the simplest addition formula for the r-function derived in 
p]], which is closely connected with the Fay's trisecant formula [21]. 

Generalized Hirota identity gives rise also to other addition formulae 
from [PJ. Indeed, let us choose in (0) 



5 (x)<rV) = n 



a=l 



v - a c 
v — b„ 



(78) 



where n is an arbitrary integer and x' — x = X)S=i[ a a] ~~ I n this case 
equation (|j) gives 



n 



M ~ a c 

\fJ--ba 



( n \ n A ci ( n \ 

X I A, n, x + J] [o a ] - n _ A, /i, x + ^ [6 a ] + 



a=l 



J - L 1 A — 6 a 



a=l 



a=l 



7 ,77^a ' 



n 



X A,6 a ,x + ^[a 5 ] x (6a,M.x + Z][ 6 *] = °-( 79 ) 



<5=1 



8=1 



Substituting the expression (f7q ) into (79) and shifting x — > x+ [A], one gets 
(fi - a a ) 



n(M-M(M-A) 



(A - a a ) 



t x+M + Qa,] r x+[A] + £> 7 ] 



7=1 



7=1 



n(A-6 a )(A-//) 



r x + M + £> 7 ] r x + [A] + X> 7 ] + 



7=1 



7=1 



, 



i (&« - M)(A - 6 a ) 6 a - 6 7 

r[x+[6 a ] + ^[a /3 ]]r(x+[A] + [ M ]+ £ [Ml =0, 



(80) 



0=1 
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It is not difficult to check (after some renotations) that equation (pC|) coin- 
cide with the Plucker's relations for universal Grassmannian manifold (see 
Theorem 3 of Q). That means, according to the Theorems 1-3 of the paper 
that r is the r function of the KP hierarchy. Note that the formulae 
(|76|), ( |80[ ) provide solutions simultaneously for the KP, mKP (dual mKP) 
and SM-KP hierarchies. 



7 Determinant formulae for transformations. 

The analytic-bilinear approach allows to represent in a compact from not 
only the integrable hierarchies but also rather general transformations acting 
in the space of solutions. 

We will consider here the transformations which are equivalent to the 
action of an arbitrary meromorphic function g(X) on the CBA function 
x(A, n, x). Let g(\) be a meromorphic function in G which has simple poles 
at the points en (i = 1,2, ...,n) and simple zeros at the points bi (i = 
1,2, ...,n), i.e. g(X) = n?=i x=fS To construct the transformed CBA func- 
tion x^(A,/i) it sufficient to find a solution of the equation 

x{v,n)g{v)xWv)dv = 

dG 

where \^ has the same normalization ((A — as x ■ 

The simplest way to find consists in the use the following consequence 
of equation (||) 

/ dvdpx{v,v)g{v)x ] {p,v,g)g~ 1 (p)x{\p) = 0. 
dG JOG 

Using this formula, one finds 

X t (A,^)= 5 - 1 (AM/,)^±i (81) 

where 

A m , a/ 3 = x( a a,b/3), a,fi = 1,2, ...,m 

and a n+ \ = A, b n+ i = fi. 

The formula (^lj) defines the generic discrete transformation of x- in 
terms of the r— function this transformation has the form 

rt(x - [A] + [ M ]) _ (A _ M)r i (AM/i) . detF n+1 



rt(x) T(x)detF n 
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where 

t(x - {a a ] + [bp]) 
a a — bp 

In the simplest case n = 1 we have (bi = b , a\ = a) 



F, 



n ,af3 



X f (A,/i) 



and 







X(a,b) x(a,/x) 


(A 


-b)(p-a) 


X(A,6) x(A,m) 


(A 


-a)(ji-b) 


X(a, b) 






r(*-[a]+ffl) 



(82) 



rt(x - [A] + [//]) (A - /x)(a - 6) (A - b){fi - a) 



.... r(x-[a] + M) 

a— b a— a 

r(x-[\] + [b]) r(g-[Al+M) 

A— 6 A— /i 



rt(x) (A-a)(/x-6) t(z)t(x - [a] + [6]) 

where = det A. One can represent the transformations ( |8l| ) also in terms 
of the function tp(X, fi) . In that form the determinant formulae ( |8l| ) taken 
at A = or /i = are very similar to the determinant formulae for the 
iterated Darboux transformations (see e.g. p2| ). The formulae (81) give us 
the Darboux transformations for all subhierarchies (KP, mKP, dual mKP, 
SM-KP hierarchies) of the generalized KP hierarchy. 

The determinant formulae ( |8ll) provide us also the multilinear relations 
for for the r-function. Since t'(x) = r(x + J2i'=i(~[ a i] + M))> the formula 
( |8lD is the n-linear relation for the T-function. It is easy to check that in the 
simplest case n = 1 this formula gives the simplest addition formula (77). 
At n = 2 the formula looks like 



r(x)r(x - [A] + [/i] + [h] - [ ai ] + [b 2 ] - [a 2 ]) x 

(A - n)(n - ai)(// - a 2 )(A - &i)(A - fc 2 ) x 
(n - bx)(n - 6 2 )(A - ai)(A - a 2 ) 

t(x-[A] + M) 

t(x+ [6i] - [a x ] + [62] - [02]) x 



r(x-[ai] + [bi]) r(x~[a 1 ] + [fc 2 ]) 

ai— 61 ai— 62 

r(x-[a 2 ] + [fei]) t(x-> 2 ] + [6 2 ]) 

0,2— b\ Cl2—b'2 



r(x-[A] + [fti]) r(x-[A] + [b 2 ]) 



A— [1 




A-61 




A-6 2 


r(x— [ai]4 


N) 


r(x-[ai] + 


[61]) 


r(x-[a 1 ] + [6 2 ]) 


ai— ix 




01— 61 




ai— fe 2 


r(x-[a 2 ]4 


M) 


r(x-[a 2 ] + 


[61]) 


r(x-[a 2 ] + [6 2 ]) 



a 2 -/t 



02— 6l 



a 2 — b 2 



(83) 



It is cumbersome but straightforward check that the relation (33) is satisfied 
due to the higher addition formulae (BO) with n = 2. 
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8 Global expression for the r-function and the closed 
one-form 



In this section we will investigate in what sense the formula ( |7q ) defines the 
r-function. Using this formula, it is possible to prove that the r-function is 
defined through the CBA function in terms of the closed 1-form, which can 
be written both for differentials of the variables x and for variations of the 
function g [g is defined on the unit circle). Moreover, it is possible to show 
that the formula ( |7G| ) and the definition in terms of the closed 1-form are 
equivalent. 

For calculation in terms of variations, it is more convenient to use the 
formula (76) in the form 

T (g(v) x 

nsW)(A -n) 

(see |jn|). Differentiating (|34|) with respect to A, one obtains 

i r * (s(y) * (&)) i , d 



dv = — (A - ix)x{\ g) , (85) 



r(g) J 5logg v-X dX 

or for A = fi 

— rx f —f] -dv = xr{X,X;g) 86 

r(g) J Slogg v-X 

where Xr{\^',g) = x(A,Mj<7) — (A — fi)^ 1 - This equation can be rewritten 
in the form 

1 fSr(g(u)) 1,1, N 1 , 

-dv=&xMv\g) rdv. (87) 



r(g) J 5\ogg v - X J ' v-X 

Relation (^7j) implies that the functionals ~^^ Sr /f ^J > and Xr{v,v;g) are 
identical for the class of functions analytic outside the unit circle and de- 
creasing at infinity. Thus for ^ belonging to this class one has 

-5\ogr(g(u)) = <f Xr (u, v; g)^-dv . (88) 

This expression defines a variational 1-form defining the r-function. It is 
easy to prove using the identity @ that this form is closed . Indeed, ac- 
cording to 

&x{\n;g) = j> x(v,^g)^^x{\v-,g)di> , 
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/S (is') 
Xr{v } \;g)— i —-Xr{\v;g)dv. (89) 

So the variation of the ( |88[) gives 

-5 2 \ogr(g) = M Xr {u, A; g) Xr (X, u; g) S ' 9 ^ ^ (A) dv dX . (90) 
// g{y) g(X) 

The symmetry of the kernel of second variation with respect to A, v implies 
that the form (|88|) is closed. 

So the formula (|88|) gives the definition of the r-function in terms of the 
closed 1-form. For the standard KP coordinates 

8g(\) _ dxn 

5(A) " h xn ' 

This formula allows us to obtain a closed 1-form in terms of dx n 



00 Qn 

8logr(g(u)) = V — Xr{v,v;g) 

n=0 



dx n , . (91) 

u=0 



For x = x\ this formula immediately gives the standard formula 

d 2 1 

^ l0gT= 2 U 

where u is a solution for the KP equation. 

In fact it is possible to prove that the function r defined as the solution 
of the relation (88) satisfies the global formula (f76|). To do this, we will show 



using the formula ( |88| ) that the derivatives of difference of logarithms of the 
l.h.s. and the r.h.s. of the expression (|84|) with respect to A , \x , A , p, and 
the variation with respect to g are equal to zero for arbitrary A , [i ,g. That 
means that l.h.s. and r.h.s. of (|76| ) could differ only by the factor, and the 
normalization of the function x implies that this factor is equal to 1. 
First we will calculate the derivative with respect to A 

d ( ( , , (v — \ 

dX 

d r ( f ... (v' — XW 1 



dX 



logx(A,M)(A-M) -logr (g(u) x {^^^j + logr(p(i/))^ 
logx(A,M)(A - fi) + j>Xr {v,v',g{v')) x {^r~j) ~\~^ dv = 

^log X (A ) M)(A-M)+Xr(A,A; ff ( i /)x (^)) (92) 
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The second term can be found in terms of x(A, ji; g) using the determinant 



formula (82). The formula fl82|) gives 



x(a,m;5 x 



v — X 



(A-//)(/i-A) 



det 



x(A,wsO x(A,//;#) 

^x(A,M5ff) |x(A,/;3) 



M — A 1 ' x(b,a;g) 
Taking the regular part of this formula at /i = A, one immediately obtains 
that the expression ( p2| ) is equal to zero. 

The case with the derivative over \i is analogous; derivatives over A and 
Jl immediately give zero. So now we proceed to the calculation of variation 

•v- A 



(93) 



8 log x(A, /i)(A - n) - log r x 



v — \i 



+ \ogr{g{v)) 



1 



x(A,/u) 



x(^ m; f)^rvx(A, i>; sO<Zf 



Xr{v,u;g(v) x )— r^r + 



x\y-,v\g)——dv 



(94) 



(we have used (|S9|) and (|88|)). Using the determinant formula (^TJ) to trans- 
form the second term, one concludes that the variation is equal to zero. 

Below we will give a brief sketch of derivation of 1-form in terms of 
Baker-Akhiezer function. Differentiation of (|76|) (with shifted arguments) 
with respect to A and \i gives 







d 



dX + 9xi (A) 



d 







dfj, dx\(n) 



log[(A - /i)x(A,^,x)] 
log[(A - /u)x(A, n, x)] 



where 
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9xi (A) 



E^iA"- 1 ^-. Equation ([ 
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logr(x) 



1 



gn-1 



implies that 



^A) logT(x) ' (95) 
(> logr(x) (96) 



777" + 







dx n '~°' y "' (n - 1)1 dX n - 1 \\dX ' dxi(X) 
Hence, we have the closed 1-form u 

lo(x , dx) = 



log[(A - /i)x(A,/i,x)] 



A=o 



E 



i 







n— 1 



\ (n-iy.dX^ 1 \\dX 9xi (A) 







+ 



d 



log [(A - m)x(A, x)] 



dx n (97) 



A=0 



and oj = dlogr. Similar expression for u> can be obtained using fl9q). At 
/i = the formula ( p7| ) is equivalent to the formula found in Q . 
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